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Abstract 

We investigate a dynamically triangulated random surface action that consists of a gaussian term 
plus the modulus of the intrinsic scalar curvature. We find that the flips are frozen out and the 
internal geometry is regularized as the coefficient of the latter term is increased. Such a term thus 
provides a way of interpolating between dynamically triangulated (ie fluid) and crystalline random 
surfaces. 



The partition function for a dynamically triangulated random surface incorporates a sum over trian- 
gulations (or more generally, polygonizations) of the surface which means that we are dealing with 
a fluid surface: 

» N-l 

^ = E/ Ud^teM-s,). (1) 

If we were motivated by investigating string theory Sg could be a simple gaussian term, 

Sg = lJ2(^'~^D'' (2) 

<y> 

arising from the discretized Polyakov string action where the X's live at the vertices of the triangula- 
tion and the sum < ij > is over all the edges. It is also possible to consider area and edge length actions 
^ in this context, but all fail to lead to a sensible continuum theory because the string tension does 
not scale Q|. 

There have also been numerous simulations and analytical treatments of fixed (triangulation) or 
crystalline random surfaces in which the sum over triangulations is dropped in favour of a single random 
or regular discretization of the surface. These can be thought of either as models in their own right of 
polymerized membranes in solid state physics or as an approximation to the dynamically triangulated 
surfaces which hopefully retain some of the salient behaviour. Various actions, many incorporating a 
"stiffness" or extrinsic curvature term |^ 

Se= ^ {l-Ui- Uj), (3) 

where rii^rij are the normals on neighboring triangles Ai,Aj, have been looked at in the solid state 
physics context. Initial simulations |^ suggested that there was a second order phase transition for 
fixed surfaces (the "crumpling transition") as the coefficient of Se was increased, though later work has 
shown that the incorporation of realistic features such as self-avoidance suppresses this ^. It is clear, 
however, that the crumpling transition does continue to exist in the absence of self-avoidance on a fixed 



surface and is second order |10 . The existence of this transition on a fixed surface has spurred many 
groups to investigate actions incorporating stiffness terms on dynamically triangulated surfaces in the 
hope of finding a viable continuum limit for a string theory The initial simulations suggested a 

second order transition, just as for a fixed surface, but more recent simulations of larger surfaces indicate 
a higher order transition at which the string tension may, however, still be scaling satisfactorily. 

Analytical work is not supportive of the existence of a transition for dynamically triangulated surfaces 
[H, suggesting that the coupling of the stiffness is asymptotically free, and this is backed up by the 
calculations for fixed surfaces where the in-plane rigidity mediates the transition. The calculations for 
the dynamically triangulated surfaces do miss out the Liouville mode and it is conceivable that this plays 
an important role in a transition if one exists, by inducing long range correlations. The status of the 
transition with alternative actions, such as the "Steiner" term suggested by Savvidy et.al. is also rather 
uncertain |p^ . 

In view of the uncertainties surrounding the behaviour of dynamically triangulated surfaces it would 
be useful to have some means of interpolating between fixed or crystalline surfaces, where the crumpling 
transition is better understood, and dynamically triangulated surfaces. A regular, flat triangulation will 
have each vertex of the triangulation surrounded by six others so a term in the action of the form 

S^J2^q,-6) (4) 
i 

where qi is the number of neighbours of vertex i, immediately springs to mind as a means of smoothing 
the internal geometry, but this is proportional to the Euler characteristic of the surface and will not 
contribute to the dynamics. Similarly, we can try 

S^J2(l^~6f (5) 
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but this is the discretization of a higher derivative Rl^ term and would not be expected to contribute in 
the continuum Umit because it is irrelevant on dimensional grounds An action of this form has been 
used to simulate the statistical mechanics and dynamics of unembedded Voronoi lattices [^6) where it 
was found that it did have a smoothing effect on the geometry, but that there was no sign of a peak in 
the associated specific heat. A third possibility is to consider 

Sfreeze^^^\qi-Q\ (6) 

i 

where the modulus sign ensures that we do not get the Euler characteristic and the normalization is 
chosen to agree with [|l3| where the term was used as an observable rather than being included in the 
dynamics. This term is non-topological and will not be suppressed by inverse powers of the cutoff, as 
equ.(^ is, because its dimension is lower. 

S freeze IS also a natural object to consider from a geometrical point of view. For a curve C embedded 
in three dimensions it was shown by Fenchel that 

- / \K\ds>2 (7) 



TT Jc 

where k is the curvature and the equality (the case of a "tight immersion") holds when C is a plane 
convex curve In we investigated a version of this for surfaces 

where the outer sum is over all the vertices of the triangulation and the inner sum is round the neighbors 
j of a node i. 4>ij is the angle subtended by the jth triangle at the ith vertex. Sught would be maximized 
for convex immersions and might therefore be expected to help smooth out surfaces, but we found that 
it was insufficient to overcome the disordering effect of area or gaussian terms. S freeze is the equivalent 
of Stight for the intrinsic geometry of the surface and in what follows we will examine the effect of this 
term on both the intrinsic and extrinsic geometry of a dynamically triangulated random surface. For 
compatibility with our earlier work we include the contribution of the measure in the action 

i 

and simulate Sg + Sm + freeze for various A. 

As we are varying the coefficient of a term in the intrinsic geometry, S freeze, we would expect most 
of the effects to be intrinsic rather than extrinsic. We therefore measure S freeze itself and the variance 
(or specific heat) term associated with it 



)2 

C = — (< Sf^^^^^ > - < S freeze > ) ■ (10) 



We also measure the mean maximum number of neighbours max((7i) and the measure term Sm to get an 
impression of the irregularity of the internal geometry as A is varied. A useful check that the simulation 
is performing as it should is provided by removing the modulus sign in equ.(^) which should then sum to 
give 47r for every configuration by the Gauss-Bonnet theorem. The extrinsic geometry is monitored by 
looking at the gyration radius X2 

X2 = 7^ r V (X^ - X^)%^g,. (11) 

m{N - 1) ^ ^ ' ^' ^ ' 

which is small in a crumpled phase and large in a smooth phase, as well as by direct inspection of 
"snapshots" of the surfaces. The expectation value of Sg can be shown to be d{N — l)/2 for all A by 
exactly the same scaling argument that is used to give this result for the gaussian plus stiffness actions, 

^On very small meshes a term of the form log(gi) which generates an B? if expanded does appear to affect the surface 
configurations if it is included with a very large coefficient, but this is almost certainly a finite size effect llS). 
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in this case because S freeze is independent of X and therefore triviaUy scale invariant. This serves as an 
additional useful check of equilibration. 

The simulation used a Monte Carlo procedure which we have described in some detail elsewhere . 
It first goes through the mesh moving the X^s, carrying out a Metropolis accept/reject at each step, and 
then goes through the mesh again carrying out the "flip" moves on the links, again applying a Metropolis 
accept/reject at each stage. The entire procedure constitutes a sweep. Due to the correlated nature of 
the data, a measurement was taken every tenth sweep and binning techniques were used to analyze the 
errors. We carried out lOK thermalization sweeps followed by 30K measurement sweeps for each data 
point for the 72 node surfaces and 50K measurement sweeps for the 144 node surfaces. The acceptance 
for the X move was monitored and the size of the shift was adjusted to maintain an acceptance of around 
50 percent. The acceptance for the flip move was also measured, but this can not be adjusted. As we 
shall see below it is this measurement that provides the greatest interest. 



A 


sweeps 


Sg 




S freeze 


C 


X2 


max{qi) 


Fs 


0.500 


30K 


106.38(0.01) 


123.63(0.00) 


99.13(0.01) 


0.37( 0.00) 


1.94(0.01) 


11.64(0.00) 


0.47 


1.000 


30K 


106.56(0.04) 


125.42(0.00) 


57.76(0.01) 


0.84( 0.00) 


2.09(0.01) 


9.27(0.00) 


0.26 


1.500 


30K 


106.48(0.06) 


126.17(0.00) 


35.08(0.04) 


1.04( 0.00) 


2.00(0.01) 


8.06(0.00) 


0.12 


1.750 


30K 


106.39(0.06) 


126.37(0.00) 


27.80(0.03) 


1.04( 0.00) 


1.88(0.02) 


7.62(0.00) 


0.08 


2.000 


30K 


106.25(0.04) 


126.52(0.00) 


22.42(0.03) 


0.98( 0.00) 


1.71(0.02) 


7.30(0.00) 


0.05 


2.250 


30K 


106.39(0.05) 


126.60(0.00) 


18.86(0.07) 


0.87( 0.01) 


1.63(0.02) 


7.09(0.00) 


0.04 


2.500 


30K 


106.62(0.05) 


126.68(0.00) 


16.21(0.01) 


0.62( 0.00) 


1.54(0.00) 


6.88(0.00) 


0.03 


3.000 


30K 


106.48(0.02) 


126.73(0.00) 


13.91(0.02) 


0.34( 0.01) 


1.54(0.00) 


6.49(0.00) 


0.01 



Table 1 
Results for iV = 72 



The results for the 72 node surfaces are shown in Table 1 . If we look at Sg we can see that it remains 
close to d/2{N — 1) = 106.5 for all the A values so the results appear to be fairly well equilibrated. We 
have not tabulated the Euler characteristic obtained by dropping the modulus sign in S freeze but this also 
comes out correctly. The results for X2 show that S freeze does not smooth out the extrinsic geometry 
as A is increased indeed, if anything, it has the reverse effect - see Fig. 1. However, the effect of the 
term on the intrinsic geometry is much more marked. We can see from the values for max{qi) that the 
geometry is regularized quite rapidly as A is increased and that by A = 3 the surfaces have an essentially 
flat internal geometry. Even more striking is that the flip acceptance Fs drops sharply to essentially zero 
as A is increased. We have not seen this with the other actions we have simulated such as Sg + XSe or 
Sg + XStight, which remain fluid at any apparent phase transitions. We have thus frozen the internal 
geometry by the addition of S freeze- 



A 


sweeps 


Sg 


Sm 


S freeze 


C 


X2 


max{qi) 


Fs 


0.500 


50K 


214.49(0.01) 


249.11(0.00) 


202.23(0.02) 


0.39( 0.00) 


2.53(0.01) 


13.01(0.00) 


0.47 


1.000 


50K 


214.42(0.03) 


252.92(0.00) 


115.57(0.04) 


0.89( 0.00) 


2.76(0.01) 


10.14(0.00) 


0.27 


1.250 


50K 


214.49(0.02) 


253.88(0.00) 


87.89(0.01) 


1.03( 0.00) 


2.93(0.02) 


9.29(0.00) 


0.19 


1.750 


50K 


214.38(0.07) 


254.92(0.00) 


52.02(0.03) 


1.10( 0.00) 


3.24(0.08) 


8.20(0.00) 


0.08 


2.000 


50K 


214.80(0.06) 


255.19(0.00) 


40.78(0.10) 


1.11( 0.01) 


3.53(0.31) 


7.83(0.01) 


0.05 


2.250 


50K 


214.30(0.04) 


255.36(0.00) 


32.20(0.12) 


1.08( 0.00) 


3.46(0.06) 


7.50(0.01) 


0.03 


2.500 


50K 


214.09(0.04) 


255.51(0.00) 


25.96(0.03) 


0.93( 0.01) 


2.85(0.08) 


7.27(0.00) 


0.02 


3.000 


50K 


214.78(0.07) 


255.69(0.00) 


17.68(0.03) 


0.62( 0.00) 


1.78(0.01) 


6.97(0.00) 


0.01 



Table 2 
Results for N = 144 



The results from the simulations of the 144 node surfaces. Table 2, are in accordance with the above. 
There is still a gentle bump in the speciflc heat associated with S freeze which gives some indication that 
we may be seeing a genuine transition rather than a crossover, though it is difficult to say for certain with 
the modest simulations we have carried out here. The value of Sg (~ 214.5 in this case) and measurements 
of the Euler characteristic again indicate that the simulation is performing as it should. Even on these 
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larger surfaces S freeze is most effective at freezing out the flip moves on the surface, which reassures us 
that we are not seeing a finite size effect. The value of max((7i) 6) indicates that the geometry is not 
completely regular in the frozen phase. The effect of the spherical topology is to reduce the average value 
of qi (an icosahedron, for example, has = 5 for all the points) so there must be defects of = 7 frozen 
in too, though there cannot be many as S freeze is close to its minimum bound of 47r. 

From these results it is clear that the addition of S freeze to a gaussian action allows us to interpolate 
between a dynamically triangulated surface with a fluctuating internal geometry and a surface with a 
fixed, regular {qi ~ 6) geometry. As this is an internal effect that appears to be largely decoupled from the 
extrinsic geometry, there is no reason to expect anything different for other actions such as those used in 
hunting for continuum string theories. An interpolation between fixed and dynamical surfaces with this 
method would then be of great interest in examining the crumpling transition as the intrinsic geometry 
is gradually frozen (or unfrozen!) to see how it changes. The intermediate state where the surface is not 
completely frozen might allow the modelling of fiuid surfaces with some rigid internal structure p| wh ich 
have previously been investigated using simulations in which a rigid frame is explicitly included pO[. It 
would also be useful to examine cases such as the Ising model where analytical results are available to 
determine the nature of the crossover from fixed to dynamical critical exponents. 

This work was supported in part by NATO collaborative research grant CRG910091. CFB is supported 
by DOE under contract DE-FG02-91ER40672 and by NSF Grand Challenge Applications Group Grant 
ASC-9217394. The computations were performed on workstations at Heriot-Watt University. We would 
like to thank R.D. Williams for help in developing initial versions of the dynamical mesh code. 



^Such as red blood cells. 
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Figure Captions 

Fig. 1. A typical surface at A = 3. There is no sign of any smoothing in the extrinsic geometry. 
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